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Finite Element Analysis of Thermal Distortion Effects on Optical Performance of 
Solar Dynamic Concentrator for Space Station Freedom 


Michael P. Doherty and Vithal Dalsania 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


Summary 

An analysis has been performed to predict the thermal 
distortion of the solar dynamic concentrator for Space Station 
Freedom in low Earth orbit and to evaluate the effects of that 
thermal distortion on concentrator on-orbit performance. The 
analysis required substructural finite element modeling of 
critical concentrator structural subsystems, structural finite 
element modeling of the concentrator, mapping of thermal 
loading onto the structural finite element model, and the 
creation of specialized postprocessors to assist in interpreting 
results. Concentrator temperature distributions and thermally 
induced displacements and slope errors and the resulting 
receiver flux distribution profiles are discussed. Results 
determined for a typical orbit indicate that concentrator facet 
rotations are less than 0.2 mrad and that the change in facet 
radius due to thermal flattening is less than 5 percent. The 
predicted power loss due to thermal distortion effects is less 
than 0.3 percent. As a consequence the thermal distortions 
of the solar dynamic concentrator in low Earth orbit will 
have a negligible effect on the flux distribution profiles within 
the receiver. 


Introduction 

To generate power for Space Station Freedom, nasa will 
use both photovoltaic and solar dynamic power-generating 
systems. Four photovoltaic modules will provide 75 kWe 
(kilowatts of electrical power) for the phase 1 Space Station 
Freedom. Twenty-five-kWe solar dynamic modules are being 
developed to provide additional power for future growth. The 
solar dynamic power-generating system is more than four times 
more efficient in converting thermal energy to electric power 
than the photovoltaic system. The closed Brayton cycle solar 
dynamic system collects solar energy to heat a working fluid, 
which in turn drives a turbine to rotate an electrical generator. 
A thermal energy storage medium is required to heat the 
working fluid during the 28- to 36-min eclipse portion of the 
95-min orbit. 

To collect and focus sunlight for a 25-kWe solar dynamic 
module, nasa will develop a large solar concentrator that 
will provide a distributed solar flux within a heat receiver 
cavity. Since reflective surfaces of this size have not been 
developed or flown in space, there is a need to predict the 
optical performance of such large, lightweight solar concen- 


trators. In order to achieve appropriate flux distribution within 
the heat receiver cavity, concentrator mirrored-surface accur- 
acy and pointing accuracy tolerances are tightly controlled. 
Identified error sources that will influence the optical 
performance of the concentrator while it is on orbit include 
facet alignment in a 1-g environment, facet manufacturing 
errors including slope error and specular reflectance, receiver- 
to-concentrator alignment, and thermal distortion. In order to 
maintain its optical performance over a 10-year lifetime, the 
concentrator must be durable enough to withstand the harsh 
environment of Space Station Freedom ’s low Earth orbit (leo). 
The leo environment causes degradation of optical and 
structural surfaces due to atomic oxygen, ultraviolet (uv) 
radiation, micrometeoroid impacts, and a considerable number 
of thermal cycles. 

This report discusses the analytical determination of solar 
dynamic concentrator thermal distortion and its effects on 
concentrator on-orbit optical performance. First, the 
concentrator is briefly described to establish a familiarity with 
its geometry. Then, the details of the analysis are presented. 
Specifically discussed herein are substructural finite element 
analyses conducted to study concentrator component structural 
behavior, the creation of a finite element structural model of 
the concentrator based on findings from the substructural 
analyses, the development of a specialized analytical procedure 
for determining thermal distortion and its effects on 
concentrator optical performance, and the execution of the 
concentrator analysis. Finally, the results of the structural 
analysis (the thermal distortion) and the effects of thermal 
distortion on concentrator on-orbit optical performance are 
presented and discussed. 

Concentrator Description 

The concentrator is one of the eight major assemblies of 
Space Station Freedom ’s solar dynamic module, as shown in 
figure 1. The concentrator is subdivided into 19 hexagonal 
panels (hexpanels) sized to fit in the space shuttle’s payload 
bay. The concentrator is an offset parabolic configuration in 
which the flat hexpanels are fastened together by latches so 
that each latch point lies on the paraboloid. The primary 
advantages of the offset parabolic and hexpanel design concepts 
are the low mass moment of inertia of the solar dynamic 
module about Space Station Freedom ’s transverse boom and 
the compactness of the stowed module, which allows for 
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packaging two complete modules in the shuttle payload bay. 
The concentrator reflective surface area comprises 456 
facets, 24 facets per hexpanel. The hexpanels are supported 
by a nine-strut support structure. Three of the struts provide 
stiffness to the hexpanels, and the other six struts attach the 
concentrator assembly to the receiver interface ring at three 
points. The receiver is tilted approximately 51° with respect 
to the boresight axis of the Sun to improve the circumferential 
flux distribution on the heat receiver cavity wall. The thermal 
design strategy for the concentrator includes using low- 
absorbance and low-emittance surface coatings to control 
component temperatures and selecting materials with a very low 
coefficient of thermal expansion to reduce thermal distortions. 

Because the solar dynamic concentrator has not yet been 
developed, the concentrator analyzed for this report was not 
the flight concentrator but rather the advanced-development 
concentrator designed and developed by the Harris 
Corporation, Government Aerospace Systems Division (ref. 1). 
nasa contracted with Harris to develop and demonstrate the 
most effective means of collecting and focusing solar energy 
to be used in a power-generating system for Space Station 
Freedom. The advanced-development concentrator assembly 
is mapped to a spherical surface rather than to the ideal 
parabolic surface. The spherical surface shown in figure 2 
allows the hexpanels to be equally spaced on a sphere. Note 
that the projected views shown in figure 2 distort the 
appearance of the equally spaced panels. The equal spacing 
of the hexpanels reduced fabrication costs by decreasing the 
number of unique latch configurations and drawings. 

The details of a hexpanel are shown in figure 3. A hexpanel 
comprises twelve rectangular-cross-section, graphite-fiber- 
reinforced epoxy (gfre) box beams. The 2-m (6.56-ft) long 
box beams are joined and bonded together at the hub and six 
corner points by shear plates and corner fitting assemblies. 


SPHERICAL MAPPING 
CONTOUR -i 



FIGURE 2. - CONFIGURATION OF ADVANCED DEVELOPMENT CONCENTRATOR. 



TYPICAL HEXAGONAL CORNERS WHERE THREE BOX BEAMS 
ARE BONDED TO CORNER FITTINGS AND SHEAR PLATES 


FIGURE 3. - DETAILS OF HEXAGONAL PANEL. 
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The gfre box beams have a high stiffness-to-weight ratio, a 
high strength-to-weight ratio, and a very low coefficient of 
thermal expansion (-0.9 x 10 6 cm/cm °C). The top and 
bottom gfre shear plates provide the load path between the 
box beams at the six comers and the hub. The aluminum comer 
fittings define the hexpanel geometry and provide the 
attachment point for the aluminum latches (ref. 2) 

A typical aluminum latch-striker assembly is shown in 
figure 4. The latch-striker assembly is a self-locking “ball and 
socket” connector comprising a striker assembly (housing and 
ball) that fits into a latch assembly (housing, spring-loaded 



pawl, and striker plate). The striker ball is locked into place 
with the spring-loaded pawl of the latch assembly. This design 
provides zero translational displacement in three axes. 

In order to concentrate the Sun’s rays into the heat receiver, 
24 mirrored triangular facets are mounted to each hexpanel 
with 72 (three each) standoff-flexure assemblies. The standoff- 
flexure assemblies isolate the facets from loads imposed by 
boxbeam distortions and allow individual facets to be aligned 
so that solar flux will be properly distributed within the heat 
receiver cavity. In order to adequately support the facets, the 
standoff-flexure assemblies are located both at the ends (near 
the corner and hub fittings) and at the midspans of the box 
beams. The details and dimensions of a typical corner-located 
standoff- flexure assembly are shown in figure 5. The facets 
have varying spherical radii of curvature and are tilted within 
the hexpanel frames depending on their specific location on 
the concentrator. The triangular facets for the advanced- 
development concentrator measure approximately 1 m on a 
side and have surface contour radii of 1921, 2181, 2441, or 
2702 cm (756.25, 858.75, 961.25, or 1063.75 in.). Four radii 
facets instead of 456 unique radii facets were chosen for the 
advanced-development concentrator in spite of a slight increase 
in slope error because of the large decrease in facet manufac- 
turing (tooling) costs. As shown in figure 6, each facet is made 
of two 7.5-mil gfre facesheets bonded to a 0.635-cm (0.25-in.) 
thick aluminum honeycomb core. The vapor-deposited reflec- 
tive and protective surfaces consist of aluminum with a 
magnesium fluoride or silicon oxide coating. Aluminum, 
although less reflective, was selected over silver for the 
advanced-development concentrator because of its durability 
in the terrestrial environment. The specular surface for the 
reflective coating is provided by an epoxy-rich layer of graphite 
vail on the facesheet and a polished-surface caul plate (ref. 3). 




FIGURE 5. - DETAIL AND DIMENSIONS OF STANDOFF-FLEXURE ASSEMBLIES. (DIMENSIONS ARE IN 
CENTIMETERS.) 
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FIGURE 6. - FACET DETAILS. 


Thermal Distortion Analysis 

The analysis to determine thermal distortion and how it 
affects concentrator performance consisted of four parts: 
substructural finite element modeling of distinctive concen- 
trator components, creation of a concentrator finite element 
model, development of a specialized analytical procedure for 
determining thermal distortion and its effects, and execu- 
tion of the concentrator analysis. The concentrator components 
modeled and analyzed as distinct substructures were a single 
triangular facet and an entire hexpanel. The findings from 
these detailed component analyses resulted in a simplified 
concentrator model. The creation of a finite element model 
of the concentrator required remodeling a single hexpanel 
with simplifications, reproducing the hexpanel at 19 separate 
locations, and adding the remaining necessary modeling detail. 
A specialized analytical procedure for determining thermal 
distortion and its effects was needed because the analysis 
had unique pre- and postprocessing requirements. Computer 
programs had to be developed at the Lewis Research Center 
to meet these requirements. These in-house programs, in 
addition to two commercially available programs, patran 
(ref. 4) and msc/nastran (ref. 5), were the basis of the 
analytical procedure for determining thermal distortion and 
its effects, (patran is a pre- and postprocessor for interactive 
creation of finite element models and graphical evaluation 
of finite element analysis results; msc/nastran is a general- 
purpose computer code for structural analysis by the finite 
element method.) Execution of the concentrator analysis 


involved running nastran and performing the postprocessing. 

The four parts of the analysis are described in detail in the 
following sections. 

Substructural Analysis 

The concentrator has two distinct components: the facet and 
the hexpanel. These two components were modeled in detail 
to study the effects of thermal loads. Analysis of these sub- 
structures was necessary to establish the degree of modeling 
detail required to accurately predict the thermal distortion of 
the concentrator. These models were created and analyzed by 
using patran and nastran. The following paragraphs describe 
the models and the results of the substructural analyses. 

Facet model . — The facet model consisted of one facet and 
three standoff-flexure assemblies. The facet chosen for sub- 
structural analysis was assumed to be supported by three 
midbeam-located standoff-flexure assemblies, a support con- 
figuration representative of one-quarter of the facets making 
up the concentrator. The other facet support configuration, 
two midbeam-located standoff-flexure assemblies and one 
corner-located assembly, was not modeled because of the 
similarity in the standoff-flexure assemblies between two 
designs. The model was used to study the effects of uniform 
temperature rise and linear through-thickness temperature 
gradient on the single facet and the standoff-flexure assembly. 
The model was also used to investigate the influence of facet 
curvature on structural stiffness. In addition, the facet was 
modeled with two different mesh densities so that the effects 
of mesh density on the accuracy of the results could be predicted. 

The finite element model of a single triangular facet is shown 
in figure 7. The composite facet structure was initially modeled 
as being flat and comprising 105 nastran tria3 plate 
elements. Seven bar elements were needed to model each of 
the midbeam-located standoff-flexure assemblies. The facet 
facesheet-honeycomb structure was modeled as a sandwich 
plate. Its membrane and bending behavior were described by 
quasi-isotropic material elastic constants (gfre facesheets), and 
its transverse shear behavior was described by isotropic 
material elastic constants (honeycomb core). Boundary 
conditions were established by assuming the points where the 
flexure assemblies attached to the hexpanel beams to be fixed. 
A temperature of 73 °C (163 °F) was applied to the standoffs 
in the form of temp cards, and the results were compared with 
cases where the standoffs were left at room temperature, 21 °C 
(70 °F). A temperature gradient of 2.2 deg C/cm (10 deg F/in.) 
from the hotter concave side to the colder convex side was 
applied to the facet with a temppI card in order to investigate 
facet-flattening behavior. Then, the facet was remodeled with 
a coarse 10-element mesh and the results were compared with 
the results obtained with fine 105-element mesh. Finally, the 
facet was remodeled to have spherical curvature, using a fine 
105-element mesh density. 

The results and conclusions of the substructural analysis of 
the facet were as follows: 
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(1) Thermal loading on the standoffs affected facet rigid- 
body rotation and therefore must be considered in the 
concentrator analysis. 

(2) The through-thickness temperature gradient on the facet 
caused an out-of-plane displacement of the center of the facet 
equal to 0.025 cm (0.01 in.) relative to the corners. This 
resulted in an increase in the radius of curvature of 58.4, 76. 1, 
96.4, and 116.8 cm (23, 30, 38, and 46 in.) for each of the 
four facet radii (1921, 2181, 2441, and 2702 cm, respectively). 
Since a change in facet radius of curvature is equivalent to 
a change in facet focal length, temperature gradients on the 
facets will have a direct effect on concentrator optical 
performance. Therefore, temperature gradients through the 
facet must be accounted for in the concentrator analysis. 

(3) Stiffening effects associated with curvature were small 
and can be neglected. Therefore, the facets in the concentrator 
finite element model can be modeled with “flat plate” finite 
elements. 

(4) Nodal displacements predicted with the coarse 
10-element mesh corresponded closely with those predicted 
with the fine 105-element mesh. Therefore, coarse structural 
modeling of the facet is adequate, and a coarse mesh can be 
used to model the facets in the concentrator model. 

Hexpanel model — The hexpanel finite element model is 
shown in figure 8. The model comprised hexpanel frame box 
beams, 24 facets with associated standoff-flexure assemblies, 
and two striker assemblies at each hexpanel corner. The 
model’s complexity reflected the degree of uncertainty in the 
structural behavior of the hexpanel. The model was used to 
study the effects of uniform temperature rise and differential 
temperature through the hexpanel. 

The hexpanel was restrained at the 12 striker ball locations 
in the out-of-plane direction only. The hexpanel frame box 


beams were modeled with quad4 plate elements so that the 
effect of differential temperature through the hexpanel from 
the “hot side” to the “cold side” could be accurately 
predicted. Also, modeling the frame beams with plate elements 
(rather than bar elements) gave depth to the hexpanel frame, 
enabling a realistic representation of the attachment of standoff- 
flexure assemblies and striker assemblies to the frame, as can 
be seen in the figure. Because facet curvature was ignored, 
each facet was modeled with four flat tria3 plate elements. 
The model comprised 804 nodes and 930 elements. Isotropic 
material behavior was assumed for all components. Two 
thermal load cases of importance were run: 

(1) Uniform temperatures of 88 °C (190 °F) on the facets 
with the remaining structure at room temperature 

(2) Differential temperature through the hexpanel frame: 
7.2 °C (45 °F) on the Sun-facing side of the hexpanel frame 
box beams and 1 .6 °C (35 °F) on the opposite side of the box 
beams, with the remaining structure at 4.4 °C (40 °F) 

The results and conclusions of the substructural analysis of 
the hexpanel were as follows: 

(1) Because the frame saw negligible movement when 
thermal loads were applied only to the facets, the standoff- 
flexure assemblies are effective in isolating the facets from 
the hexpanel frame. The reason for this desired behavior is 
the relative flexibility of the standoff-flexure assembly with 
respect to the hexpanel frame. Therefore, it is possible to 
greatly simplify the modeling details required to represent the 
standoff-flexure assembly in the concentrator finite element 
model. It is also possible now to accept the increases in facet 
radius of curvature due to thermal gradients predicted by the 
facet substructural model as being indicative of and applicable 
to all the facets in the concentrator. As a result, temperature 
gradients on the facets will not be revisited in the concentrator 
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STANDOFF-FLEXURE ASSEMBLY 
FOR THREE ADJACENT FACETS 

FIGURE 8. - HEXPANEL FINITE ELEMENT MODEL. 


finite element analysis, but the known increases in facet radii 
of curvature due to thermal gradients will be accounted for 
directly in the nasa Lewis in-house programs that predict 
thermal distortion effects. 

(2) Differential temperature through the depth of the hex- 
panel frame box beams caused large out-of-plane displacements 
to the hexpanel, indicating that temperature differentials (or 
gradients) through the hexpanel frame should be considered 
in the concentrator analysis. Because bar elements can model 
temperature differentials across both the depth and the width, 
they can be used to model the hexpanel frame box beams in 
the concentrator finite element model. 

The conclusions drawn from the substructural analysis 
enabled the finite element representation of the concentrator 
to be simplified. The modeling of the concentrator based on 
the findings from the substructural analyses is discussed in the 
next section. 

Concentrator Finite Element Modeling 

The finite element model of the solar concentrator was 
created by using patran. All modeling information was stored 
in a patran data base that could be updated and revised as 
required before being translated to nastran bulk data form. 


A single hexpanel was remodeled with simplifications based 
on the findings from the substructural analyses. The remodeled 
hexpanel appears in figure 9. The hexpanel contains 223 grid 
points, 204 bar elements, and 72 triangular plate elements. 
Each hexpanel frame box beam was modeled by using four 
equal-length bar elements attached end to end. The bar 
elements were connected in all six degrees of freedom at the 
hexpanel corner and center node points. Four bar elements 
per box beam was fine enough modeling fidelity to allow 
accurate attachment of standoff-flexure assemblies at the 
corners and midsections of the frame, because of the simplified 
representation of the standoff-flexure assemblies in the 
concentrator model. 

Each standoff-flexure assembly now became represented by 
a flexible bar element, which, in order to have its proper 
location with respect to the hexpanel frame maintained, was 
attached to the hexpanel frame corner or midsection by a rigid 
bar element. (This geometry can be seen in the blowup view 
in figure 9, which shows one continuous box beam, six rigid 
bar elements, and five flexible bar elements (the sixth is 
hidden) having unique local coordinate systems.) The flexible 
bar element was given a stiffness equivalent to that of the entire 
physical standoff-flexure assembly. The derivation of this 
equivalent stiffness is shown in the appendix. 
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FLEXIBLE BAR 
ELEMENTS HAVING 
EQUIVALENT STIFFNESS 
OF STANDOFF-FLEXURE 
ASSEMBLY 


FIGURE 9. - CONCENTRATOR FINITE ELEMENT MODEL - SINGLE HEXPANEL. 


The facets were modeled coarsely and without curvature by 
using flat triangular plate elements. Three plate elements repre- 
sented each facet, enabling a node to be located at the facet 
center to give a more direct readout of facet rotations and 
translations. In the actual design a standoff is attached to the 
facet approximately 7.0 cm (2.75 in.) from the facet corner. 
So that the standoffs could be connected at the facet corners 
in the concentrator finite element model, the actual facet 
corners were shifted toward the center of the facet to the 
approximate standoff- facet connection point. Consequently, 
the facet corners were located only approximately in the finite 
element model and had to be later modified with the actual 
locations for postprocessing by the nasa Lewis in-house 
programs. 

The concentrator model was constructed by using patran 
commands to reproduce this hexpanel at 19 separate locations. 
The 19 panels had to be interconnected to form a single 
structure. The panels were interconnected by adding latch and 
striker details at the proper locations. Figure 10 shows the 
concentrator finite element model including the latch and 
striker modeling details. Each latch-striker assembly was 


modeled by using four flexible bar elements and four rigid 
bar elements. The junction point between the latch and the 
striker, being a ball joint, was represented by releasing the 
rotational degrees of freedom at the four bar ends connected 
at the junction grid point. The rigid bar elements were 
necessary to provide offsets between the panel corner grid 
point and the latch-striker bar ends because the hexpanel frame, 
being of bar element representation, was without cross- 
sectional depth. 

The concentrator is supported by nine tubular struts. Three 
struts are used to provide in-plane rigidity for the concentrator; 
the remaining six struts are used to attach the concentrator 
to Space Station Freedom at the solar receiver. Bar elements 
were used to model the tubular struts. The junctures of the 
struts and the concentrator are pinned connections. The three 
nodes at the solar receiver interface are fixed in all six 
directions. 

The patran modeling of the concentrator was then 
completed by adding the material property cards and the bar 
element property cards (tables I and II). The material and the 
optical performance requirements of the concentrator 
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TYPICAL LATCH-STRIKER ASSEMBLY 


RIGID BAR 
ELEMENTS 

FLEXIBLE BAR 
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RIGID BAR / 
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FIGURE 10. - CONCENTRATOR FINITE ELEMENT MODEL INCLUDING LATCH-STRIKER DETAILS. 


TABLE I.— MATERIAL PROPERTIES 


Material location 

Modulus of 
elasticity, 
MPa 

Poisson’s 

ratio 

Coefficient of 
thermal expansion, 
cm/cm °C 

Aluminum standoff- 
flexures and latch- 
strikers 

6.895 x 10 4 

0.25 

3.600X10- 6 

Composite hexpanel 
box beams 

1.793 xlO 5 

.30 

-.900x10-6 

Honeycomb facets 

3.650X10 3 

.30 

3.600x10-6 

Graphite tubular 
struts 

2.275x105 

.30 

.018x10-6 


TABLE II.— BAR ELEMENT PROPERTIES 


Component 

Cross-sectional 

area, 

cm 2 

Moment of 
inertia, 

l r 

cm 4 

Moment of 
inertia, 

Iv 

cm 4 

Torsional 

constant, 

cm 4 

Hexpanel box 
beam 

2.6128 

40.4519 

3.1150 

9.6803 

Standoff- 

flexure 

.1826 

2.418XI0- 3 

5.2x10-5 


Latch-striker 

1.9613 

8.7409 

.4162 

1 .4152 XlO 2 

Tubular strut 

9.5768 

76.807 

76.807 

1.5359X10 2 


demanded a more effective characterization of the effects of 
thermal distortion than could be inferred directly from the 
patran nodal displacements. The additional finite element load 
case preparation tasks, as well as the need to predict and 
display the effects of “thermally induced slope error’’ on 
concentrator optical performance were undertaken and 
addressed by the nasa Lewis in-house computer programs. 
These in-house computer programs are described in the 
following section. 

Analytical Procedure 

Four in-house computer programs were used to perform 
finite element load-case preparation tasks as well as to evaluate 
results from the finite element analysis. These four computer 
programs were panel, snip (ref. 6), conrms, and offset 
(ref. 7). panel and snip were used to complete the 
preparations for the finite element load case: panel revised 
the locations of the facet corners within the concentrator model 
to accurately reflect the Harris design coordinates; snip 
mapped temperature results defined under a parallel thermal 
analysis effort (ref. 8) to the finite element model, conrms 
and offset were used to evaluate results from the finite 
element analysis: conrms made an approximate optical error 
calculation for each facet; offset performed a ray-tracing 
analysis to predict actual changes in flux distribution at the 
solar receiver due to the thermally distorted concentrator 
geometry. Figure 1 1 shows a flow chart of the entire analytical 
procedure, with the program name displayed in bold, the input 
and output of each program shown, and the in-house programs 
marked with an asterisk. The following subsections describe 
these in-house computer programs in further detail. 
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NAME OF PROGRAM 


INPUT/OUTPUT 


PANEL 


FACET CORNER COORDINATES 
FROM DRAWINGS 


PATRAN/PRE 


PATRAN/PATNAS 


SINDA-NASTRAN INTERFACE 
PROGRAM (SNIP) 


NASTRAN 


FACET GRID AND ELEMENT CARDS 


NASTRAN BULK DATA CARDS 
REPRESENTING CONCENTRATOR 
FINITE ELEMENT MODEL (FEM) 


TEMPERATURE, GRADIENTS, 
GEOMETRIC DATA, AND NODE 
NUMBERS FROM THERMAL ANALYSIS 


TEMPERATURE LOAD CARDS 
FOR NASTRAN ELEMENTS 


DISPLACEMENT VECTORS 


PATRAN/NASPAT 

L 

L— , 1 

r 


REFORMATTED DISPLACEMENT 
VECTORS AND ROTATIONS 
FOR PATRAN PLOTS 





ADDITIONS OF RMS COLUMN 
TO PATRAN PLOT FILE 


DISPLACED COORDINATES OF 
FACET CORNERS 


OFFSET |t< 

1 




FLUX DENSITY ON RECEIVER 
WALLS AND APERTURE PLANE 


PATRAN/POST USING FOUR 
DIFFERENT MODELS 


r * 
i 

-I 

I 

L 


TEMPERATURES AND NODE NUMBERS I 
FOR PATRAN CONTOUR PLOTS FROM I 
THERMAL ANALYSIS 1 

1 


CONTOUR PLOTS OF DISPLACEMENTS, 
RMS ERROR, TEMPERATURE 
DISTRIBUTIONS, AND RECEIVER 
FLUX DISTRIBUTION 


FIGURE 11. - CONCENTRATOR THERMAL DISTORTION ANALYSIS PROCEDURE. (ASTERISK DENOTES PROGRAMS 
DEVELOPED AT LEWIS RESEARCH CENTER.) 
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PANEL.— The. corner locations of the facets in the finite 
element model were updated by using the Fortran program 
panel, panel accepts facet*, y, and z comer coordinates from 
the Harris drawings and produces new nastran connectivity 
cards for the facets. This file of connectivity cards is then 
translated into neutral file format and merged with the original 
patran data base containing the concentrator model. This 
merging procedure causes the original facet node and element 
information to be overwritten with new facet node and element 
information and the corner locations of all 456 facets to be 
revised to correctly reflect the Harris design coordinates. 

panel was necessary because the concentrator modeling 
process, being based on the reproduction of a generic hexpanel, 
had created facets having corner locations that were only 
approximate in the depth direction within the hexpanel frame. 
Having accurate facet corner locations in the finite element 
model was important for two reasons. First, accurate location 
of the facet corners imposes a precise length to each standoff 
(flexible bar element) in the finite element model. A precise 
representation of standoff physical length was necessary 
because substructural analysis had indicated that thermal 
loading on the standoffs must be considered in the concentrator 
analysis. Second, the offset code required displaced facet 
corner locations in order to predict flux profiles for distorted 
concentrator geometry. The process of obtaining displaced 
facet corner locations is simplified by initially having accurate 
facet corner locations in the finite element model. 

SNIP.— The sinda to nastran Interface Program (snip) was 
used to eliminate the labor-intensive task of relating 798 
specific thermal nodes of a parallel thermal analysis to the 4400 
specific nastran structural elements of the concentrator 
model. The thermal analysis had been performed by using the 
Thermal Radiation Analysis System (trasys) (ref. 9) and the 
System Improved Numerical Differencing Analyzer (sinda) 
(ref. 10). snip is a Fortran computer program that generates 
nastran thermal load cards for nastran plate, shell, bar, 
and beam elements. Inputs to the snip program are sinda (or 
similar thermal analyzer) nodal temperature results, their 
corresponding nodal coordinates, and the nastran structural 
model in bulk data file form, snip uses a geometric search 
routine and a numerical coding scheme to relate thermal model 
nodes to structural elements, snip then calculates element 
temperatures on the basis of the weighted average of thermal 
node temperatures related to each element. User-controlled 
input parameters provide control over node-to-element 
correlation. 

The concentrator was assumed to be in a 250-n mi altitude 
circular orbit, which has a period of 93.7 min with 35.7 and 
58.0 min of eclipse and sunlight, respectively. The trasys 
and sinda programs were used to predict the temperatures of 
the concentrator at moments in time during this orbit. Three 
load cases were considered in the structural analysis, corre- 
sponding to three points in time during orbit: 6, 30, and 54 
min into the sunlit portion of the orbit (approximately morning, 
noon, and dusk). For each load case the input to the snip 


program consisted of 798 node data cards, representing 24 
facets, 12 box beams, and 6 corner fitting assemblies for each 
of the 19 hexpanels. Associated with each node was a 
temperature and the centroidal x, y, and z coordinates for the 
component (facet, box beam, or corner fitting assembly) 
it represented. The node data cards representing the box 
beam cards contained a temperature equal to the average 
temperatures of the four sides of the box beam and temperature 
gradients in both the beam depth and width directions based 
on the temperature differences between top and bottom and 
left and right box beam sides. Each box beam had a local 
coordinate system in order to achieve the correct orientation 
of the temperature gradients within the structural element. 

Certain types of thermal loading were identified from 
substructural analysis as being important to consider in the 
concentrator structural analysis: temperatures for the facets, 
the box beams, and the standoffs; and temperature gradients 
for the facets and the box beams. After user selection of the 
appropriate qualification and numerical coding parameters that 
relate the thermal and structural models, snip generated 
nastran element thermal load cards for the facets, the box 
beams, and the standoffs for each load case. For each standoff 
bar element snip produced a weighted average temperature 
based on the available temperatures of the nearest facet and 
hexpanel corner fitting assembly. For each box beam snip 
identified with each of the four bar elements making up that 
box beam a temperature and both a depthwise and a widthwise 
temperature gradient. For each facet snip related a uniform 
temperature to all three plate elements making up the facet. 
Temperature gradients on the facets were accounted for 
directly in the ray-tracing analysis. 

CONRMS. — A Fortran program called conrms was written 
to produce root-mean-square (rms) values of the rotation vector 
for all 456 facets. This rms rotation, an approximate optical 
error calculation, was required because the nastran output 
displacement vector is useful only for describing gross 
concentrator deflection and does not give any sure indication 
of individual facet misalignment due to thermal distortion. 

For each facet conrms calculates an rms rotation by squaring 
the nastran displacement vector rotations about the * and 
y axes (the bending rotations) at all four nodes of the facet, 
summing these values, taking the square root of the sum, and 
dividing the sum by four. These rms rotations are written on 
the sixth column of the nastran displacement vector for the 
respective grid points. These rms rotations are then able to 
be displayed graphically with patran alongside the nodal 
displacements, as a more sophisticated indicator of optical 
performance degradation due to thermal distortion. In like 
manner, conrms calculates rms rotations for each panel and 
for the entire concentrator for all three load cases. 

conrms also calculates the displaced facet corner locations, 
needed as input into the offset ray-tracing program. Displaced 
locations are calculated by adding the *, y, and z deflections 
of the facet corner grid points to the actual coordinates of each 
of the 456 facets as supplied on the Harris drawings. 


10 


OFFSET . — A ray-tracing computer code called offset was 
developed at the nasa Lewis Research Center specifically to 
predict incident flux profiles for the offset solar collector of 
the Space Station Freedom solar dynamic electric power 
system. This program traces rays from 50 points on the face 
of the Sun to 10 points on each of the 456 collector facets and 
from the facets through the solar receiver aperture to the walls 
and back face of the receiver. The facet corner locations are 
the principal input to the program. If the design facet corner 
locations are input, the resulting solar flux distribution on the 
receiver surface characterizes the undeformed orientation of 
the facets. If displaced facet corner locations due to thermal 
distortion are input, the effect of concentrator thermal 
distortions, or thermally induced slope error, on receiver flux 
distribution can be determined. The offset program also 
predicts the effect of the two other optical error sources: 
changes in facet radius of curvature due to thermal flattening, 
and statistically distributed (rms) surface slope error. 

Graphical output of offset results requires a patran model 
of the receiver that uses plate elements to display the energy 
per unit area on each of the receiver walls. This patran model 
appears in figure 12. The receiver is cylindrical with a small 
aperture centered on one end to admit reflected solar rays. 
It has a diameter of 1.86 m (6.10 ft) and a length of 2.99 m 
(9.81 ft). The cylindrical surface was divided into 80 elements 
in the tangential direction and 18 elements in the axial 


direction. The back face was divided into 80 circumferential 
by 6 radial elements; the front aperture plate was divided into 
80 circumferential by 5 radial elements. The offset code 
calculates the flux density for each plate element and creates 
the patran file necessary to plot the receiver flux distribution 
maps for the deformed concentrator. The offset code also 
calculates the power lost to the receiver by summing the power 
on the elements making up the front aperture plate. 

Execution of Analysis 

Once the modeling revisions of panel and the thermal load 
preparation by snip were completed, the three nastran 
structural load cases to determine concentrator thermal 
distortion were run. Figure 13 shows the thermal loading on 
the concentrator facets for load case 2 at 30 min into the sunlit 
portion of the 250-n mi altitude orbit (noon). Although not shown 
here, thermal loading simultaneously existed on the hexpanel 
box beams and the standoffs for this (and each) load case. 

Once the nastran results were obtained, conrms was run 
to calculate the rms rotations and the displaced facet corner 
locations, offset used the displaced facet corner locations to 
calculate the receiver flux distribution profiles. As a final step, 
all graphical results (temperatures, displacements, rms 
rotations, and receiver flux distribution profiles) were 
displayed on patran. 



FIGURE 12. - PATRAN MODEL OF SOLAR RECEIVER. 
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FIGURE 13. - THERMAL LOADING ON FACETS FOR LOAD CASE 2 (NOON). 


Discussion of Results 

A contour plot of total displacements of the concentrator 
facets for load case 2 is shown in figure 14. Note that total 
displacements varied widely across the structure. Figure 15 
displays, in a head-on view of the concentrator, the facet rms 
rotations for the same load case. The results produced by 
conrms indicate that rms rotations for the entire concentrator 
are 0.04, 0.12, and 0.17 milliradian (mrad) for load case 1 
(morning), load case 2 (noon), and load case 3 (dusk), 
respectively. Even the maximum concentrator rms rotation of 
0.17 mrad is small, only about 6 percent of the 3-mrad facet 
slope error fabrication tolerance. This indicates from a 
preliminary standpoint that the effects of thermal distortion 
on concentrator on-orbit optical performance are small. 

Results from the offset ray-tracing analysis are presented 
in table III. This table shows the power lost to the receiver 
for each thermal load case as well as that lost from facet surface 
slope error and facet thermal flattening. It can be seen that 
a fabrication tolerance of 3 mrad for the facets (as is expected 
practically) was responsible for the greatest amount of spillage. 
A comparison of rms calculations and optical analysis results 
shows that rms calculations appear to be a good qualifying 
indicator of the effects of thermal distortion on concentrator 
performance. The combined case of 3-mrad slope error, 
thermal flattening of the facets, and thermal distortions actually 
decreased the amount of power lost from the 3-mrad error 
case. The reason is that the largest facet radius of curvature, 
2702 cm, was actually undersized by about 200 cm (78 in.). 


As the facets flattened from thermal distortion, the increased 
radius of curvature became closer to the ideal radius of 
curvature. In other words, the image at the aperture was made 
smaller rather than larger. 

Figure 16 shows the flux distribution profile on the receiver 
walls for the combined case. No significantly high flux rates 
were observed, with the maximum being 3.5 W/cm 2 . The 
flux distribution maps were identical for all practical purposes 
for the three thermal load cases and the perfect facet case. 


TABLE III.— POWER LOST DUE TO DIFFERENT 
OPTICAL ERRORS 


Case description 

Power lost 

kW 

Percent of 
total 

Perfect facets (undeformed orientation) 

0.4685 

0.226 

Facets with 3-mrad slope error 

8.3713 

4.042 

Change in radius due to thermal 

.3081 

.149 

flattening of facets (0.0254 cm) 



Thermal load case 1 (6 min into Sun) 

.4515 

.225 

Thermal load case 2 (noon) 

.4730 

.228 

Thermal load case 3 (3 min before 

.4770 

.228 

eclipse) 



Combined thermal load case 2, thermal 



flattening, and 3-mrad slope error 

7.9308 

3.830 
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FIGURE 15. - ROOT-MEAN-SQUARE ROTATIONS OF CONCENTRATOR FACETS FOR LOAD 
CASE 2. 
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FIGURE 16. - RECEIVER FLUX DISTRIBUTION PROFILE FOR COMBINED CASE - ON-ORBIT 
THERMAL LOADS WITH AN INITIAL 3-mrad FACET SLOPE ERROR. 


Concluding Remarks 

A finite-element-based analysis has been performed on the 
Space Station Freedom solar dynamic concentrator to predict 
its performance in the 250-n mi low-Earth-orbit thermal 
environment. A finite element structural model of the 
concentrator has been developed. On-orbit temperature 
distributions from the thermal analysis have been used for 
thermal loading. Specialized optical analyses, making use of 
the finite element displacement results, have assisted in result 


interpretation. The results show that the distortions of the 
concentrator are extremely small with respect to the size of 
the structure and that the effect of thermal distortions on 
concentrator optical performance should be negligible. In fact, 
the effect of on-orbit concentrator thermal distortions on solar 
dynamic system performance was found to be less than that 
caused by fabrication tolerances for optical facet hardware. 
This indicates the choice of materials for the solar dynamic 
concentrator to be excellent from the thermal design 
standpoint. 
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Appendix— Derivation of Equivalent Stiffness for Standoff-Flexure Assembly 


Each standoff-flexure assembly in the concentrator finite 
element model is represented by a flexible bar element attached 
to the hexpanel frame box beam by a rigid bar element. This 
flexible bar element has a stiffness equivalent to that of the 
entire standoff-flexure assembly. The derivation of this 
equivalent stiffness follows. 

Figure 5 shows the details and dimensions of the corner- 
located standoff-flexure assembly. The reference coordinate 
system is shown in this figure. As shown 

L length of flexure bracket, 4.303 cm (1.693 in.) 

T thickness of flexure bracket, 0.076 cm (0.030 in.) 

W width of flexure bracket, 0.953 cm (0.375 in.) 

H height of standoff, 5.080 cm (2.000 in.) 

D diameter of standoff, 0.483 cm (0.190 in.). 

Let 

def y = As y + Af y + (H0 fz ) 
and 

def z = As z + Af z + (H6 fy ) 

where 


P load applied at standoff ball in y direction, N 
H height of standoff, 5.080 cm 

E Young’s modulus for aluminum, 6.895 XlO 4 MPa 
(10.0X 10 6 psi) 

I sz moment of inertia for standoff about z-z axis, 
*D 4 /6 4, 2.6715 XlO -3 cm 4 (0.653 XlO -4 in. 4 ) 

Therefore, 

As y = 2.372 x 10 “ 3 P 


From reference 11, page 97, case Id: 

PL 3 


4f,= 


192 £/ a 


where 


(A2) 


P load applied at flexure midspan in y direction, N 
L length of flexure bracket, 4.303 cm 
E Young’s modulus for aluminum, 6.895 XlO 4 MPa 

Ifx moment of inertia for flexure about x-x axis, UT 3 /12, 
3.49X10 -5 cm 4 (0.838X10 -6 in. 4 ) 


def y total deflection at standoff ball in y direction due 
to load P 

def z total deflection at standoff ball in z direction due 
to load P 

As y deflection of standoff ball with respect to standoff 
base (jam nut) in y direction due to load P at 
standoff ball in y direction, cm 

As z deflection of standoff ball with respect to standoff 
base (jam nut) in z direction due to load P at 
standoff ball in z direction, cm 

Afy deflection of flexure bracket in y direction due to 
load P at its midspan in y direction, cm 

Af. deflection of flexure bracket in z direction due to 
load P at its midspan in z direction, cm 

dfy rotation of flexure bracket about y axis, rad 

Ofi rotation of flexure bracket about z axis, rad 

Calculation of Def y 


Therefore, 

Afy = 1.724 XlO -3 P 
From reference 11, page 287, eq. (1): 

* t( 

e * = J5 (A3) 

where 


t twisting moment, (P/2)//, N-cm 

2 length of member, or half-length of flexure bracket, 
L/2, 2.152 cm (0.847 in.) 

J polar moment of inertia of flexure bracket, 

(W72)(772) 3 [(16/3) -3.36(77110], 1.324xl(T 4 cm 4 
(0.318X 10 -5 in. 4 ) 

G modulus of rigidity for aluminum, 1.793 XlO 4 MPa 
(2.6 XlO 6 psi) 

Therefore, 


def y = As y + Afy + (H6 fi ) 


6^ = 2.3025 XlO -2 P 


From reference 11, page 96, case la: 


ASy = 


PH 3 

3 K 


As a result 


(Al) 


def y = Asy + Af y + (H6 fi ) 

= [2.372 + 1.724 + (5.080x23.025)] XlO -3 P 


where 


= 0.121063 P 
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Calculation of Def z 

def z = As z + Af z + (HOjy) 

Since the standoff bar is of circular cross section, 

As z = As y = 2.372 x 10 -3 P 

The flexure bracket is assumed to have no deflection in the z 
direction owing to a load P at its midspan in the z direction; 
therefore, A f z = 0. 

From reference 11, page 103, case 3d, the following 
expression can be derived: 


\6EI fy 

where 

M applied couple, PH = 5.08 P, N-cm 
L length of flexure bracket, 4.303 cm 
E Young’s modulus for aluminum, 6.895 Xl0 4 MPa 
I fy moment of inertia for flexure about y-y axis, 7TT 3 /12, 
5.4816X 10 -3 cm 4 

Therefore, 


only calculate the principal flexural moments of inertia for 
a 5.080-cm-long cantilevered bar that yield the same y 
and z deflections (def y and def z) calculated here. From 
equation (Al): 


def y = 


pi if 

3EL eq 


where 


H = 5.080 cm 
E = 6.895 X10 4 MPa 


def y = 0.121063 P 


As a result 

heq 


^(5.080) 
3£(0. 121063 P) 


5.2x 10~ 5 cm 4 


Likewise, 

_ P(5.080) 3 
Iyeq ~ 3 E def z 

where 


E — 6.895 Xl0 4 MPa 


- 3.6X10" 5 P 


def z = 0.002555 P 


As a result 


As a result 


def z = As. + Af z + (HOfy) 


I yeq = 2.481 X 10' " 3 cm 4 


= [237.2 + 0.0 + (5.080x3.6)] xl0~ 5 P 
= 0.002555 P 

To make the flexible bar element have a stiffness equivalent to 
that of the entire standoff-flexure assembly, one need 


In the concentrator finite element model a local coordinate 
system is defined for each flexible bar element so that these 
principal flexural moments of inertia, I yeq and I zeq , are of 
proper orientation. In the blowup view in figure 9, local 
coordinate systems for two of the six flexible bar elements 
are shown. 
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